Electron transport through asymmetric ferroelectric tunnel junctions: 

current-voltage characteristics 
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We have carried out calculations of current-voltage characteristics for the electron tunnel current 
through a junction with a thin insulating ferroelectric barrier assuming that interface transmissions 
for the left and right interfaces noticeably differ due to dissimilarity of the interfaces. Obtained con- 
ductance vs voltage and current vs voltage curves exhibit well distinguishable asymmetric hysteresis. 
We show that the asymmetry in the hysteretic effects could originate from the asymmetric bias volt- 
age profile inside the junction. In particular, we analyze the hysteresis asymmetries occurring when 
the bias voltage distribution is low-sensitive to the spontaneous polarization reversal. 

PACS numbers: 73.23.-b, 85.50.-n 
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Ferroelectric materials attract significant interest of 
the research community due to their potential useful- 
ness in various technological applications [ij. In the 
past decade it was repeatedly shown in the experiments 
that ferroelectricity could be maintained in films of some 
perovskite oxides with thickness of the order of a few 
nanometers These experimental observations were 

corroborated by first-principles computations predicting 
the critical film thickness for ferroelectricity to persist in 
perovskite films to be as thin as a few lattice parame- 
ters [3-[l3l- Thus, both experimental and theoretical re- 
sults indicate that ultrathin ferroelectric films exist and 
may be used as barriers in ferroelectric tunnel junctions 
(FT J). A ferroelectric tunnel junction consists of two 
conducting electrodes (leads) separated by a nanometer- 
thick ferroelectric film. Electrons can tunnel across the 
film thus maintaining a nonzero conductance of the junc- 
tion. Specific nature of the ferroelectric barriers add new 
functional properties to FTJs (Tlj Such junctions could 
serve as essential components in manifold nanodevices 
such as binary data storage memories and switches [l^ . 

The most important feature of ferroelectrics is the pres- 
ence of spontaneous polarization which could be reversed 
by an applied electric field. The polarization switch- 
ing inherent to ferroelectric materials alters the sign of 
polarization charges at the electrode/barrier interfaces 
thus changing the electrostatic potential profile inside the 
junction. Also, the polarization switching shifts ions in 
ferroelectric unit cells and changes lattice strains in the 
film affecting the electron band structure of the junction. 
As a result, transport characteristics of a FT J may be 
significantly modified. The effect of the polarization re- 
versal on FTJ conduction was predicted in earlier works 
(see e.g. Refs. flsUlB]) and observed in experiments 
on Pt/Pb{Zro.52TioAs)03/SrRu03 FTJ [11, on junc- 
tions with Lao iBiQ CjMnO^ tunnel barrier [l6(] and on 
Pt/SrTiOs/Pt junctions 

Currently, theoretical and computational efforts are 



mostly concentrated on studies of zero bias conductance 
of FTJs. However, the theoretical analysis of current- 
bias voltage (/ — V) relationships is also important 
for it provides more thorough characterization of the 
electron transport through FTJs. As shown in Ref. 
[l3| . the polarization reversal in the ferroelectric films 
causes a hysteretic behavior of both current-voltage and 
conductance-voltage characteristics. Especially interest- 
ing hysteretic effects may occur in asymmetric FTJs 
where the interfaces are dissimilar. The dissimilarity 
could occur when the leads are made out of different met- 
als 14, 1^. Also, it could originate from the atomic and 
electronic structure of the junction [3]. When interfaces 
differ, the effects of the polarization reversal also may 
differ depending on the specific properties of the inter- 
faces. This may result in asymmetry of tunnel electrore- 
sistance (TER) versus bias voltage curves and of hystere- 
sis loops in current-voltage characteristics. In strongly 
asymmetric FTJs both TER and hysteresis could be well 
pronounced for a certain bias voltage polarity and sig- 
nificantly reduced on the polarity reversal. The asym- 
metric TER and hysteretic effects were recently reported 
for asymmetric junctions including BaTiO^ 2^, 21 1 and 
Pb{ZrQ,2Tia.s)03 [2§] ferroelectric barriers. 

In the present work we theoretically analyze transport 
characteristics of an asymmetric FTJ. We simulate the 
FTJ by a ferroelectric film of thickness d sandwiched 
between two semi-infinite metal electrodes. We assume 
that the film is uniformly polarized in the direction per- 
pendicular to the electrode planes, which is the direction 
of "z" axis within the chosen coordinate system. At 
small values of the applied bias voltage the conductance 
of a tunnel junction per area A could be described by 
the following expression [23| : 
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where k| | is the projection of the tunneling electron wave 



2 



vector to the electrode's plane, and E is the tunnel en- 
ergy, which takes on values close to the Fermi energy Ep. 
The transmission function T{E, k||) may be factorized 
and presented in the form 24]: 

T(i;,k||) = tL(i;,k||)exp{-2K(i;,k||)d}t,j(£;,k||) (2) 

Here, ti^{E,'k\\) and tpi{E,'k\\) are the interface trans- 
mission functions for the left and right electrodes, re- 
spectively, and k||) is determined by the height and 
shape of the potential barrier in the tunnel junction and 
by the energy-momentum relation of a tunneling elec- 
tron. Generally speaking, it also depends on the bias 
voltage but one may disregard this at small values of V. 
Employing a simple form for the energy-momentum re- 
lation: 
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where m is the effective mass of the tunneling electron, 
we may approximate: 
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Here, ^{z) is the overall potential profile seen by the 
tunneling electrons. At zero bias voltage the potential 
<i>(z) is a superposition of the potential, which deter- 
mines the positions of conduction bands in metal elec- 
trodes with respect to the Fermi energy, and the bar- 
rier potential occurring due to the presence of the insu- 
lating ferroelectric film. Also, the overall potential in- 
cludes a contribution created by surface charges on the 
surfaces of the ferroelectric film and screening charges 
appearing on the electrode surfaces. The simplest ap- 
proximation for '^{z) appropriate for a symmetric FTJ 
is a trapezoid barrier with a thickness d and heights 
$(0) = U + (y5(0), = U + ip{d) where ip{z) is the 

potential induced by the surface charges. The relative 
heights of the barrier at z = and z = d are deter- 
mined by the potential difference (p{0) — f{d). The dif- 
ference changes its sign when the polarization in the film 
reverses causing the hysteresis in the conduction through 
the FTJ to occur. In asymmetric FTJs the change in the 
complex band structure induced by the left-right asym- 
metry of the ferroelectric displacements in the barrier 
may noticeably change the average barrier heights U for 
two polarization states, as was theoretically shown for 
a SrRuOs/ BaTiO^/ SrRuOs junction basing on first 
principles computations [3]. The difference in the bar- 
rier heights for different polarization directions may be 
of the same order or greater than that induced by the 
potential (p{z). In the following transport calculations 
we take into account both the effect of depolarizing field 
corresponding to the potential fiz) and variations in 
the barrier height appearing due to the junction asym- 
metry. However, to avoid extra complications in further 



computations we do not include in our model the effect 
of an intrinsic electric field throughout the barrier, which 
appears due to the difference in the electrostatic dipoles 
at the interfaces, as well as the strain effects due to piezo- 
electricity. 

To estimate the current through the junction we use 
the expression: 



J 
A 



2e 



J dEiU-fn) J ^r(i?,k||) 



(5) 



Here, the transmission function r(i?, k||) is given by Eq. 
(m and /l.r are Fermi distribution functions with chem- 
ical potentials ij-l,r- The chemical potentials are shifted 
with respect to Ep when a nonzero bias voltage V is 
applied across the junction, namely: 

fiL^Ep + iqeV; HB, = Ep - {I - i])eV. (6) 

The division parameter 77 shows how the voltage V is 
distributed in the junction. The expression ([5]) resem- 
bles the well known Landauer formula for tunnel current 
through molecular junctions and /or quantum dots cou- 
pled to the conducting leads [25|]. It could be derived 
using surface Green's functions formalism presented in 
the earlier work [2^. We have grounds to belive that 77 
takes on values depending on the relation of the interface 
transmissions II and tn. Calculating the current den- 
sity we hypothesize that 77 = tji/{tL +tR), which resem- 
bles the expression for the division parameter commonly 
used in the theory of electron tunnel transport through 
molecular bridges and quantum dots. In that theory 77 is 
described by the similar expression where the parameters 
determining the coupling of a molecule/quantum dot to 
the left and right leads take on the parts of the interface 
transmissions. 

Regardless of symmetry /asymmetry of a FTJ, the 
transmission coefficients and Ir differ, and the dif- 
ference in their values may be significant. For instance, 
it was reported that in a Pt/BaTiO^/Pt junction one 
of the transmission coefficients was approximately three 
times greater than another one [lo| . One may expect 
even greater difference in the interface transmission for 
the FTJ studied by Maksymovich et al [l^l, where recti- 
fying current-voltage characteristics were reported. One 
reason for the disparity in the transmission functions is 
the trapezoid-like profile of the potential barrier. The 
barrier height decreases along the polarization direction, 
therefore Ir should exceed when the polarization 
points to the right and vice versa. In symmetrical FTJs 
with identical interfaces one may expect the transmis- 
sion coefficients simply interchange their values on the 
polarization reversal: 
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where an arrow indicates the polarization direction. In 
such a case the bias voltage division parameter 77 in the 
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Eq. is to be replaced by 1 — when the polarization 
direction is reversed. Accordingly, the current-voltage 
curves associated with different polarization states (al- 
though asymmetric by themselves due to the asymmetry 
in the bias voltage distribution across the junction) are 
arranged in the V — I plane in such a way that hysteretic 
features remain symmetrical. The TER versus voltage 
curve must be symmetrical upon bias voltage reversal, 
as well. 

In asymmetric FTJs the relationship between the 
transmissions and tji for different polarization direc- 
tions is more complicated. Due to an asymmetric defor- 
mation of the ferroelectric potential profile the relations 
^ do not hold any more. Besides, dissimilarity in the in- 
terfaces may affect the transmissions due to the reasons 
unrelated to ferroelectric properties of the film linking 
the electrodes. Even within the simple one-dimensional 
model of a FTJ employed in the present work, one must 
take into account specific features of the contact between 
the electrodes and the film, which could significantly dif- 
fer on the left anf right sides of the junction. For instance, 
in the case when the electrodes are made out of differ- 
ent metals, this may significantly influence the interface 
transmissions. Also, FTJs studied in the recent works 
[20l . consisted of an epitaxial ferroelectric film grown 
on the top of conducting electrode and a sharp metal tip 
serving as another electrode. We may expect a nonepi- 
taxial contact between the tip and the film to notably 
modify the corresponding interface transmission for both 
polarization directions. Some other effects, which could 
modify interface transmissions in asymmetric FTJs (and, 
in consequence, the bias voltage distribution across the 
junctions) are discussed elsewhere (see e.g. Ref. 
and references therein). So, the bias voltage distribu- 
tion across an asymmetric FTJ for different polarization 
orientations could hardly be predicted basing on some 
general considerations. For each particular FTJ the rela- 
tion of the interface transmissions (which greatly affects 
the bias voltage distribution) must be separately estab- 
lished, and the results may significantly vary depending 
on the FTJ characteristics. For instance, we remark that 
in the asymmetric FTJ whose transport properties were 
studied by Gruverman et al [2^ , < tj^ for both polar- 
ization directions. This follows from the reported shape 
of the I — V characteristics. 

In summary, various asymmetric features in the hys- 
teretic behavior of the tunnel current and conductance 
in asymmetrical FTJs originate from asymmetries in the 
tunnel barrier profile for different polarization directions 
and from specifics of the bias voltage distribution (deter- 
mined by the relationship between the interface trans- 
missions). The relative effects of these two factors could 
vary but neither one should be disregarded on general 
grounds. Here, we concentrate on the analysis of trans- 
port properties of a strongly asymmetric FTJ where the 
bias voltage distribution exhibits a relatively low sensibil- 
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FIG. 1: (Color online) Conductance ^ = ^ of the model 
asymmetric FTJ vs bias voltage for the paraelectric (dash- 
dotted line) and ferroelectric (solid and dashed lines) states 
of the film. The curves are plotted using Eqs. (Il])-(|6]) at T = 
30K, d = 1.8nm, U = O.GeV (low resistance ferroelectric 
state), U — 0.8eV (high resistance ferroelectric state), for 
II ~ StR (left panel) and in — 3tL (right panel). 



ity to the polarization reversal. Accordingly, we accept 
that II /tji remains the same regardless of the polariza- 
tion direction. 

Computing the conductance and current for our model 
FTJ we used the results reported in Ref. for 
Pt/BaTiO^/Pt junctions. So, we assumed that d = 
l.Snm, and zero bias conductances per unit cell area 
take on values G/A « 17.0 x 10~^e^/h (for the paraelec- 
tric state of the film) and G/A « 2.9 x 10~'^e'^/h (for 
its low resistance ferroelectric state), respectively [2'n\. 
We did carry out calculations for tfl/t^ — 3 {rj — 0.75) 
and tL/tji = 3 (77 0.25) for ferroelectric states of 
the film. For the film in the paraelectric state we as- 
sumed = tij. Using these values we estimated the 
interface transmissions for paraelectric and ferroelectric 
barriers in the junction. As shown in the Fig. 1, the 
conductance through the junction with the ferroelectric 
barrier should exhibit rather well distinguishable hystere- 
sis. The hysteresis is asymmetric with respect to V = 0, 
and we observe the significant asymmetry in the TER. 
At tL = 3tfl the TER is small at positive bias volt- 
age but it significantly increases upon the reversal of the 
bias voltage polarity, as shown in the left panel of the 
figure. When tn = St^ the TER is much better pro- 
nounced at positive bias voltage (see the right panel of 
the Fig. 1). Basing on the experimental results of Gar- 
cia et al [2l|] for thin BaTiO^ films, we may roughly 
estimate the coercive voltage Vc for our model FTJ as: 
Vc 1.0- 1.5 V^. Therefore, the difference in the TER val- 
ues at V = ±Vc, which occurs due to the specifics of the 
bias voltage profile across the junction may reach values 
of the order of 10^ — 10'^. The current- voltage charac- 
teristics are presented in the Fig. 2. We see asymmetric 
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FIG. 2: (Color online) Current-voltage characteristics for the 
model ferroelectric tunnel junction. The curves are plotted 
using Eqs. H])-® The parameters take on the same values 
as in the Fig. 1. 



hysteretic behavior of the tunnel current. As well as for 
conductance, the asymmetry in the hysteresis loops to a 
considerable degree originates from the bias voltage pro- 
file in the junction. We remark that the I ^ V curves 
in the right panel of the Fig. 2 resemble those recently 
reported in the work (20| . 

In conclusion, we have demonstrated that current- 
voltage and conductance-voltage characteristics of an 
asymmetric FTJ may exhibit a pronounced hysteretic 
behavior of a special kind whose distinctive feature is 
a well pronounced asymmetry in the hysteresis loops for 
different bias voltage polarities. The asymmetry in the 
hysteresis largely originates from the asymmetric bias 
voltage division in the junction. In strongly asymme- 
tryc junctions the dissimilarity of the interfaces may af- 
fect interface transmissions so much that the effects of 
the polarization reversal may be surpassed. In this case, 
the relationship between the transmissions is low sen- 
sitive to the polarization revesal, and the bias voltage 
distribution across the junction rather weakly depends 
on the polarization direction. Our results give grounds 
to expect that the bias voltage reversal in the consid- 
ered strongly asymmetric FTJs should bring significant 
changes to their transport characteristics, and this could 
be used in designing of nanodevices. 
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